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Abstract 

In an R-parity violating supersymmetric scenario, the lightest neutralino \\ is no longer a 
stable particle. We calculate the branching ratio for the decay mode Xi — * v l which occurs 
at the one-loop level. Taking into account bilinear as well as trilinear lepton number violating 
interactions as the sources of R-parity violation, we make a detailed scan of the parameter space, 
both with and without gaugino mass unification and including the constraints on the neutrino 
sector from the recent Superkamiokande results. This study enables one to suggest interesting 
experimental signals distinguishing between the two types of R-parity breaking, and also to 
ascertain whether such radiative decays can give rise to collider signals of the type 77 + /E from 
pair-produced neutralinos. 



PACS NOS. : 12.60.Jv, 13.10.+q, 14.80.Ly 



1 Introduction 



Decay of a heavy particle with a photon in the final state can often lead to experimental signals for 
physics beyond the Standard Model. Recent literature contains many studies, both experimental 
and theoretical, devoted to signals of the type 77 + ft, 7 + ft etc. |Q, in the context of various 
new physics scenarios. A rather large fraction of such studies are concerned with supersymmet- 
ric (SUSY) || theories. For example, in models with gauge mediated supersymmetry breaking 
(GMSB)0, ||, the decay of the lightest neutralino (x?) (next-to-lightest supersymmetric particle 
in these models) into a photon and a gravitino can give rise to 77 + /E (or 7 + / E) which can have 
substantial event rates even after subtracting the Standard Model backgrounds. The signal 77 + 
/E is in fact the most probable discovery channel of GMSB. Side by side, radiative decay of the 
second lightest neutralino (x®)) m a scenario based on N = 1 supergravity (SUGRA), into a photon 
and the lightest neutralino (x?, the lightest supersymmetric particle in this class of models) can 
sometimes give rise to similar final states in colliders from a pair of x^'s or fr° m a X2X1 pair. It 
requires a careful investigation to distinguish the latter type of signals from the former to establish 
the reality of GMSB vis-a-vis the SUGRA-type mode of supersymmetry breaking ||. 

In view of this, it is important to know whether there can be other new physics options which 
allow radiative decay of a heavier particle into a photon and another invisible particle, leading to 
signals of the same kind as those discussed above, and, in some cases, faking the GMSB signals 
which are being looked for so carefully. It will be obvious that the closest kinematical resemblance 
to the GMSB signal is borne when the invisible particle is a near-massless neutral fermion, for which 
the immediate candidate is a neutrino. Now, radiative decay of a heavier neutral fermion into a 
neutrino becomes a distinct possiblity in R-parity violating SUSY |7| where the lightest neutralino 
can decay in this vein at the one-loop level. This motivates one to undertake an exhaustive study 
of the radiative decay Xi — > u l i n an R-parity violating scenario, to see whether in some region 
of the parameter space this kind of a decay can have an appreciable branching ratio in spite of 
the allowed tree-level decay channels of the lightest neutralino. In this paper, we have attempted 
such a study, taking into account R-parity violation through both bilinear and trilinear terms in 
the superpotential. 

It should be mentioned that Xi — * v l has been calculated in some earlier works || where 
either the composition of the lightest neutralino has been confined to certain limits, or some specific 
regions of the SUSY parameter space have been adhered to. Our purpose, on the other hand, is 
to make a detailed scan of the parameter space, and the results presented by us highlight those 
particular regions where the branching ratio can be of the largest magnitude. This, we feel, is 
necessary to establish the authenticity of, for example, GMSB signals during collider searches. We 
also go beyond the assumption of gaugino mass unification at the scale of grand unified theories 
(GUT) and present some results with the SU(2) and U(l) gaugino masses treated as independent 
parameters, to probe whether radiative neutralino decays with enhanced rates ensue upon freeing 
the SUSY standard model from a GUT embedding. And lastly, we take into account the fact that 



R-parity violating SUSY is being invoked in recent times to explain the generation of neutrino 
masses || in consonance with the Superkamiokande data on atmospheric neutrinos 0]. Thus a 
section of our results pertains to that particular region of the parameter space where the indicated 
hierarchy of neutrino masses, together with large angle mixing between the second and the third 
generations, is reproduced. 

The paper is organised as follows. In section 2, we describe the basic framework of the R-parity 
violating model used in our calculations. In section 3, some details of the loop calculation are 
presented. Section 4 contains our numerical results for various scenarios as well as for different 
choices of the parameters. We summarise and conclude in section 5 . Detailed forms of the loop 
integrals are outlined in the appendix. 

2 Basic Framework 

The fact that all the particles in the standard model carrying baryon (B) and lepton (L) number are 
fermions, with fixed gauge multiplet stuctures, implies that B and/or L cannot be violated by one 
unit. Thus one cannot write down renormalizable terms in the Lagrangian which violate B or L or 
both. This is no longer true in the case of super symmetric theories where the particle spectrum is 
now doubled and baryon number and lepton number are assigned to the supermultiplets. One can 
now have scalar particles containing B or L, and it is possible to have terms with AL or AB = 1. 
However, the scenario has to be made consistent with the non-observation of proton decay which 
requires the simultaneous violation of B and L. This is ensured in an over-restrictive manner in the 
minimal supersymmetric standard model (MSSM), by imposing a discrete multiplicative symmetry 
called R-parity defined as R = (_i) L + 3B + 25 ) which equals +1 for Standard Model particles and 
-1 for the superpartners. An immediate consequence of R-parity conservation is that the lightest 
supersymmetric particle (LSP) is stable. 

The conservation of R-parity, however, is not prompted by any strong theoretical reason, and 
theories where R is violated through nonconservation of either B or L are perfectly consistent with 
stability of the proton. Such scenarios can be studied by genaralising the MSSM superpotential to 
the following form : 

W = Wmssm + Wft (1) 

with 

Wmssm = Ah 2 + ti^L^E] + h%Qill x b] + h%QiH 2 Uf (2) 

and 

Wft = \ijkL; L LjEl + X' ijk LiQjDl + \'l jk U^D C jD c k + €iLiH 2 (3) 

In this paper we shall assume that the B-violating piece A" is not there. Also, we discuss the 
L-violating terms in two separate categories for the convenience of analysis, considering, in turn, 



Wft with either the bilinear (eiLiH 2 ) |Tl]) 12, [b|] or the trilinear (A- and A') |l4| terms existing in 
the superpotential at a time. 
Case (1). W n = e 3 L 3 H 2 

We simplify the analysis here by assuming a bilinear R-parity violating term involving only 
the third leptonic generation. All the imporatnt features of such a scenario can be seen within 
simplified framework. 

At first sight it appears that the e-term can be rotated away from the superpotential by properly 
redefining the L 3 and H\ fields, so that 

• ' ^TT (4) 



£, = -e 3 Hi + fiL 3 
Now as a result of this rotation the superpotential takes the form: 



A 1 

- ^ftL^-^U^ + ^L^ (6) 

fJL [M fl 

which means that the A- and A' type terms are generated in general and in the last two terms 
i = 1,2. 

But more importantly, it should be noted that the scalar potential has terms of the following 
forms in the original basis: 

Vscal = Ml 3 L 2 3 + m\R\ + B x[ iH x H 2 + B 2 eL 3 H 2 + ^s^i + (7) 

where we have written down the F-term and soft breaking terms relevant for our purpose here. 
It is evident that such terms will generate, in general, a vacuum expectation value (vev) v 3 for 
the sneutrino. The basis rotation regenerates such terms, and therefore the sneutrino vev is in 
general non- vanishing as the special consequence of bilinear R-parity violation whenever soft SUSY 
breaking is there p^, 16]. 



Thus there are two ways of parametrizing a bilinear R-parity violating scenario: (i) with e-term 
present in the superpotential together with a vev of the sneutrino, and (ii) with e-term rotated 
away from the superpotential but the vev of the sneutrino (in the rotated basis) embodying the 
L-violating effects that the former would imply in, say, the neutralino and chargino mass matrices. 
Of course, the parameter e takes refuge in this case in the scalar potential together with the soft 
breaking parameter B 2 . B 2 can be eliminated by using the conditions for electroweak symmetry 
breaking, if e and v 3 are used as independent variables. Our calculations here are done in a basis 
where both e and i>3 (vev of the sneutrino) are present in the superpotential. 



Now, the presence of e and v 3 will induce a kind of mixing in the fermionic as well as in the 
scalar sector of the theory, which is typical of bilinear R-parity violation. In the fermionic sector 
neutralinos will mix with the tau-neutrino and the charginos, with the tau lepton. Consequently, 
the neutralino mass matrix takes the following form: 
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where the successive rows and columns correspond to (H2, H%, —1W3, —iB, v T ). Here 



sin/3 (cos/3) 



M and M' are the SU(2) and U(l) gaugino mass parameters respectively, /U, the Higgsino mass 
parameter, and g = \j g 2 + g' 2 ■ Evidently, this will result in the generation of (Majorana) masses 
for the third generation neutrino at tree-level through a see-saw type mechanism [[H]]. 
In such a situation, the (3x3) chargino mass matrix is 
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-gvi 
-9V3 




where v\ = (Hi), 



v 2 



(H 2 ), v 3 = (9 T ), tan/? = % and / 
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Sl, M being the 



SU(2) gaugino mass parameter. Here we have assigned (— iW~ , ,t~l~) along the rows and 
(— iW + , H2 , t~r + ) along the columns. 

Similarly, in the scalar sector neutral and charged scalar mass matrices are enlarged as a result 
of the mixing between charged sleptons and charged Higgs and sneutrino with the neutral Higgs. 
In the neutral scalar sector the scalar mass-squared matrix in the basis {Re(Hi), Re(H2), Re(v T )} 
takes the form: 



/ ml + 2Ac + 4\vf 



Mi 



-A\v\V2 + B\\i 



(10) 



4AuiU3 + /ie 

-4Au \V2 + B\[i m\ - 2Ac + 4Au| -4\vzv 2 + B 2 e 
V 4Xviv 3 + fie -4Xv 3 V2 + B2€ m 2 ~ T + 2Ac + AXv 2 . 
In a similar way, the neutral pseudoscalar mass-squared matrix becomes 3 x 3 in the basis 
{Im(Hi), Im(H2), Im(v r )} and takes the form: 
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In the charged scalar sector the charged scalar mass-squared matrix becomes 4 x 4 in the basis 
{H 2 ,H 1 ,f L ,f R } and is given by 

/ r - \g ,2 c -B ± n + \g 2 v x v 2 -B 2 e + ^g 2 v 2 v 3 -efvi \ 

2 _ -Bin + \g 2 viv 2 a + \g' 2 c + fv\ fie + ^^Vua - lf 2 viv 3 -efv 2 - Afv 3 

-B 2 e + \g 2 v 2 v 3 fie + \g 2 v x v 3 - \fv x v 3 p + \gH + \g' 2 c fifv 2 + 

V -efvi -efv 2 - Afv 3 fifv 2 + Afvi q - \g l2 c + f 2 v 3 2 J 

(12) 

with 

2 , 2/ 2 , 2 i 2\ 

r = m 2 + -5 («! +v 2 +v 3 ) 

2 , 1 2/ 2 i 2 2\ 
2 i f 2 2 

p = + j v 1 

2 i j-2 2 

q = m n + f v i 

t = {-v 2 x + vj + vi) 

I 2 2 i 2 \ 

A = (5 2 + 5 ,2 )/8 

In addition to constraints coming from electroweak symmetry breaking, the requirement of 
electric charge conservation subjects the parameters in the potential to appropriate conditions in 
any given scenario. 

In the calculations here, we make the following choice of independent parameters: 
{//, tan/3, e, v 3 , m 2 L = m? fl = = m§, A t , A),, A T , B{\. We also take the ramaining scalar 
mass parameters to be as follows: 

2 2 2 2 2 

mf- = m~- = mr- = m?~ = m~, 
t\ t2 bj b 2 i 

and 

m l,A = m l 

where, tan/3 is the ratio of two Higgs vacuum expectation values, B\ is the usual B-term 
corresponding to the /i-term, A's are the trilinear soft breaking terms. 

The bilinear R-parity violating parameters e and ^3 can be constrained from different experimen- 
tal considerations such as the tau-mass measured within the existing errors, and, most importantly, 
from the constraints on neutrino masses. The latter can be subjected to current laboratory bounds, 
but the restriction becomes more stringent on using the Superkamiokande results on atmospheric 
neutrinos. In this light, the values of e and ^3 have to be constrained in such a way that the 
sneutrino vev is less than a few hundred KeV's in the basis where e is rotated away. [] 

1 Of course, bilinear terms LiHi with i = 2, 3 have to be included in order to explain — v T oscillation. This 

requires a straightforward extension of the formalism described above. For more details, the reader is referred to Ref. 
0- 



The above kind of mixing in the chargino, neutralino and scalar sectors results in physical states 
which are superpositions of neutral(charged) leptons and neutralinos(charginos) on one hand, and 
Higgses and sleptons(sneutrinos) on the other. This implies that the Yukawa and gaugino coupling 
terms, written in terms of the physical states, will give rise to all the interactions of the A and A'- 
types. In addition, the fact that gaugino couplings play a role here makes it possible to free some 
of the terms from constraints due to gauge invariance of the superpotential. One such interaction 
term has particular importance in our calculation, namely one involving a top quark, a stop and 
a state which is dominantly a neutrino. This is clearly a result of neutrino-neutralino mixing and 
the quark-squark neutralino interaction in the Lagrangian. For more details the reader is referred 
to @. 

Also, a scenario like this allows one to have tau-neutralino-W and neutrino-neutralino-Z cou- 
plings. This results in additional decay modes of a neutralino which are characteristic of bilinear 
R-parity violation fl8|| . 

Case (2). W n = \ ijk UL 3 E% + W ijk UQpi 

In such a case the requirement of gauge invariance implies that A^s are antisymmetric in the 
first two indices due to SU(2) invariance but A^- fc s are not subjected to such constraints. Thus 
we have nine A and twenty seven A'-type couplings other than the MSSM parameter space. These 
trilinear couplings can also be constrained in various ways like, for example, lepton universality 



violation [19], neutrinoless double beta decay [20], majorana mass of neutrino |21{| , flavor changing 
neutral current processes [22] etc. Their implications in various high energy collider experiments 
such as LEP p3|, HERA [p|| or the Fermilab Tevatron have also been widely explored. In 
addition, they can be responsible for neutrino masses generated at the one-loop level pG|], 

Before we end this section, it may be remarked that although constraints on the R-violating 
parameters are frequently talked about, in practice these constraints will always have to stipulate 
a given set of values of the R-conserving (MSSM) parameters. 



3 The one-loop calculations 

The one-loop calculation of the decay Xi(p) ~~ * v(ki) + 7(^2) has been performed in the nonlinear 
R-gauge. A similar calculation for the radiative decay of heavier neutralinos into lighter ones can be 
found in Ref. [27]. In this gauge, choice of the gauge fixing term modifies certain types of vertices 
compared to, say, the 't Hooft-Feynman gauge. For example, the W + G~ r f vertex (where is the 
charged Goldstone boson) is abesnt in the non-linear R-gauge, whereas the W + W~j vertex gets 
correspondingly modified. Also a set of diagrams involving a Z — 7 transition (with an off-shell Z), 
which potentially contributes to the process under investigation, gives zero contribution when one 
sums over all the loops. 

The diagrams which ultimately contribute in the bilinear R-parity violating scenario are the 
triangle diagrams shown in Fig. 1. A subset of these will be present when one considers only 



trilinear lepton number violating terms in the superpotential, the relevant diagrams being (a) and 
(b) involving only the down-type quark/squark or lepton/sleptons in the loop. This is actually a 
consequence of SU(2) invariance of the superpotential and can be seen explicitly if one expands 
the A- and A'-type terms into their component forms. The appearance of additional diagrams in 
the bilinear R-parity violating case can be attributed to mixing in the scalar and lepton-chargino- 
neutralino sectors, once electroweak symmetry breaking takes place (see section 2). In fact, as has 
already been mentioned, one of the main reasons for undertaking this study with bilinear R-parity 
violation is the presence of the additional diagrams involving the top quark and the W-boson, which 
are potential sources of enhancement of the decay amplitude. 

Since the photon is on-shell, U(1)em gauge invariance demands that the transition amplitude 
be proportional to cx^A^e* . Thus one expects the following form for the matrix element: 

M = ig^uik^iPR - r] u 7] 1 P L )a^k2^lu(p) (13) 

where M^o is the mass of Xi, Vi and f]v are the signs of the mass eigenvalues of x± and the neutrino 
respectively. Thus, depending on r\\ and 7] u , the effective x\ v l interaction is either proportional to 
c^ 1 '^2^% or Tscr^A^e*. The radiative decay width of Xi is then given by 

g 2 ~ M% 

r(xS - »l) = (14) 

where we have neglected the mass of the neutrino and <?|o„ 7 > containing details of the loop integrals, 
has the dimension of inverse mass-squared. 

In order to calculate <?^W one has to evaluate the triangle graphs. For a given set of internal 
particles there are two diagrams which differ from each other by the direction of charge flow in 
the loop. In the figures we have shown only one set of diagrams. One must be careful about the 
diagrams involving the reverse flow of charge in the loop because in those cases one encounters 
vertices with clashing arrows. 

The integrals which appear during the loop calculations are regularized using dimensional regu- 
larization. Since there is no tree-level coupling of the type x\ v l the divergences must cancel among 
each gauge invariant subset of diagrams. For example, the diagrams (g) and (h) and their partner 
graphs where the direction of charge flow is reversed form a gauge invariant subset and it can be 
shown that the infinities cancel within this set of diagrams. Similarly, one gets a finite result out of 
the loop diagrams involving fermion/sfermion loops, i.e., the diagrams {(a), (&)} and their partner 
graphs. This is true, pairwise, for the loops {(c), (d)} and {(e), (/)} involving charged Higgs boson 
and charginos and the charged Goldstone bosons and the charginos. 

The decay amplitude shown in Eqn.(13) can be decomposed in the following way (with only 
the appropriate contributions retained, as stated earlier, with trilinear R- violating terms); 



M=Mt+M 2 (15) 
The matrix element due to the four diagrams (a), (b) and their partner graphs is given by 



M 1 = -jL-^uih^PR-wxP^fafuip) 

x J2 e f C A( A L B R ~ ArBl^M^I 1 - I 3 )] + m f (A L B L - A R B R )I 2 } (16) 
/ 

where Qf is the charge of the fermion / in units of e (e > 0), and Cf is the color factor for the 
particles in the loop. One can use this result to include contributions from all possible fermion/scalar 
loops. Thus contributions from charged Higgs and Goldstone bosons are also present in the sum 
over /. We have also explicitly used the sum over the three possible chargino states (xt) i R 
the bilinear R-parity violating scenario. One must be careful with the masses and couplings as 
well as about the electric charges of the fermions in the loop. The integrals P are expressible in 
terms of the general expressions for one- loop three-point functions |2q] , in forms that are presented 
in the appendix, where we also give the expressions for the combinations {AlBl — ArBr) and 
(AlBr — ArBl) for all the possible cases. 

In the case of stop and sbottom we have included the effect of left-right mixing. Also, the 
calculation includes possible absorptive parts of the loop integrals, which can be present in particular 
when the decaying neutralino is heavier than the W. 

In a similar way, the contributions to the matrix element from diagrams (g) and (h) and their 
partner graphs are given by 



M 2 = —^uik^iPR-wxP^fafuip) 

x J2i( A L B L ~ ArBr^M^I 1 - J 2 -I 3 )] + 2M k (A L B R -A R B L )J 2 } (17) 
k 

In both the above expressions, the quantities Al, Ar, Bl and Br are all assumed to be real. 

The loops involving the top quark (Fig. 1(a) and Fig. 1(b)) , which are present only in the case of 
bilinear R-parity violation, always seem to have a rather important effect 0. Matching contributions 
also come over a large area of the parameter space from loops involving the W-boson (and the tau 
and/or the lighter chargino, depending on the relevant mixing angles) in the propagator (diagrams 
(g) and (h)). The very presence of these diagrams causes the rates for radiative decays to be larger 
in cases with bilinear R-violating effects than in those with only trilinears, except in those where 
there is destructive interference between the two types of graphs. 

2 Our detailed scan of the parameter space, however, also reveals regions where lighter fermions contribute com- 
parably because of favoured mixing angles, much in the same way as the charm quark makes dominant contributions 
to the real part of the box diagrams for A" — K° mixing. 



4 Numerical Results 



The calculation of the branching ratio for Xi ~^ v l involves the computation of decay widths of 
the lightest neutralino in all the relevant tree-level two-and three-body decay modes. The different 
types of mixing that have been discussed in section 2 for bilinear R-parity violation open up the 
following two-body decay modes, as and when kinematically allowed: 

Xi — ► t^W , Xi —>■ v T Z, Xi — > v T h where h is the lightest Higgs boson. 

In addition, there are three body tree-level decay channels with both bilinears and trilinears. 
While the channels x? ~~ > e i"^j^fc;Xi — > v idjdk,x^i ~~ *■ e i v j e h are available in both the cases [29], 
additional decays such as Xi ~^ u iUiV and Xi ~~ vvv can take place in the former, and can be 
important when two-body decays are not kinematically allowed. 

In all the numerical results presented here, we have used At = A^ = WGeV . Our finding is that 
the branching ratio of the radiative decay is rather insensitive to the values of these parameters. 

Let us first concentrate on bilinear R-parity violating scenario. We will consider two different 
cases, namely: (i) when the supersymmetry breaking gaugino mass parameters are unified at the 
Grand Unification scale, and (ii) when SU(2) and U(l) gaugino mass parameters denoted by M2 
and Mi respectively are treated as unrelated. 

Due to reasons explained below, significant contributions to the radiative decay can come only 
when bilinear R-parity violating effects are present. There again, the relevant parameters e and the 
sneutrino vev(s) can be constrained under two different considerations. If the neutrino-neutralino 
mixing process in such a scenario is the main mechanism for the generation of (Majorana) neutrino 
masses for explaining the Superkamiokande (SK) results on atmospheric muon neutrinos, then the 



parameters have to be constrained in such a way that the quantity v' = J v fi + v r ^ s ^ ess than about 
100 KeV in a basis where the e-parameters are rotated away from the superpotential. The numerical 
results modulo such constraints are denoted by "with SK" in the corresponding figure captions, 
and the results correspond to both the second-and third-generation neutrinos in the final state of 
the radiative decays. On the other hand, if one ignores the SK constraints, then, with just the 
laboratory bounds on the tau neutrino mass, the values of e and v% (i.e. the R-violating parameters 
in the third generation) can be as large as on the order of GeV's. We have also presented such 
results, considering only the third (tau) neutrino to be there as the decay product. 

Figures 2-9 (10) show the results of a scan of the parameter space carried out in the different 
scenarios mentioned above, calculated with bilinear (trilinear) R-parity violation. The first thing 
that we note is that in all the results presented, a specific set of masses for the squarks and sleptons 
are assumed. The branching ratios of the radiative decay, however, are rather insensitive to their 
variation. This is because if, for example, we reduce the squark masses, the effect of a lighter stop 
will increase the decay width for the radiative decay, but the corresponding tree level decays will 
also undergo a boost via diagrams mediated by squarks. 

Similarly, a comparison of figures 2-4 with 5-7 and 8 with 9 reveals that imposing the SK 
constraints does not cause any appreciable reduction to the probability of obtaining branching 



ratios on the higher side for similar combinations of MSSM parameters. All that it does is to 
scale the overall coefficients instrumental in both the loop-and tree-level decays. This results in 
branching ratios of similar magnitudes, but is manifested in larger decay lengths for the neutralino 
when R-parity violating parameters are restricted to yield Am| ~ a few times 1CP 3 . 

As has already been mentioned, the major contributions to the radiative decay come from (a) 
the top-stop loop with bilinear R-pariy violation, as well as from loops involving the W-boson and 
the tau or the lighter chargino. However, this leads in some regions of the parameter space to the 
interesting possibility of their cancelling each other. When such a cancellation is very severe for 
some specific neutralino mass, the branching ratio is seen to undergo a sharp dip at that point, as 
seen in Fig. 3 and Fig. 8. In such cases, the effective contributions at those points hardly contain any 
input that is special of the bilinear terms. On the other hand, in Fig. 2 (and partially in Fig. 3) one 
encounters a situation where the contribution is mainly from the W-loops but the latter undergoes 
a destructive interference among the various component terms, causing the overall branching ratio 
to fall at a particular region. 

The numerical results show a sensitivity to the MSSM parameters \x and, to a somewhat lesser 
extent, on tan (3. It is also clear from the graphs that the radiative decay tends to remain suppressed 
for small values of which is treated here as a free parameter. It gradually rises with and 
then almost saturates, showing a slight fall for approaching a TeV. The loops have the largest 
contributions for ~ 500GeV. By and large, the corresponding regions of the parameter space 
have the lightest neutralino almost entirely dominated by the Bino state. 

When the condition of gagino mass unification is relaxed, the mass parameters Mi and Mi, 
corresponding to the U(l) and SU(2) gauginos can be unrelated [30|. Most of the features of the 
unified scenario, including the possibilities of having branching ratios close to 10 per cent, are seen 
here also. However, it entails the additional possibility of having a destructive interference between 
the top-and W-induced diagrams in cases where the branching ratio is otherwise on the higher 
side, causing the latter to fall sharply by about 5 orders of magnitudes for a particular value of the 
neutralino mass, with all other parameters at same values. This is seen when M2/M1 is smaller 
than what it would have been with the constraint of unification, and the dip is found to occur when 
the two masses are quite close to each other. 

Figure 10 is a sample showing the order of magnitude of the radiative decay with only trilinear 
R-violating couplings present in the theory. The values of the A- and A'-type coupling constants have 
been used consistently with the current limits |3l]]. We have already seen that the two potentially 
most important classes of diagrams, namely, those mediated by the top and the W, are absent in 
such a case. Therefore, it is hardly surprising that the branching ratios cannot rise above 10~ 8 - 
10~ 10 over a large part of the parameter space. The range in which the branching ratio tends to 
lie in such a case matches, as one would expect, with the one to which it falls when a cancellation 
of the large contributions takes place in the scenario with bilinears. In any case, the radiative 
branching ratio turns out to be too small to be of any observable significance where R-parity is 



violated only through trilinear terms. 

With the maximum value of the branching ratio for the radiative decay being between 5 and 
10 per cent, two-photon signals from such decays will fake signals like those of GMSB in only a 
rather small region of the parameter space. Such a thing might happen when the lighest neutralino 
in GMSB is close to the kinematic limit of production. However, the fact that only the bilinear 
R-violating terms can boost the branching ratio to the level of close to 10 per cent has rather 
interesting implications in terms of observing distinctive signals of the latter. For example, pair- 
produced neutralinos at LEP energies (assuming an intergated luminosity of 500 pb^ 1 ) can give 
rise to about 40 events where there is a radiative decay on one side, leading to signals of the type 
3/ + 7. There can be many more of such events in a high-energy electron-positron collider with an 
integrated luminosity of about 50 With a proper event selection strategy, radiative neutralino 

decay can thus be an interesting signature of R-parity violation with bilinear terms. 

5 Summary and conclusions 

We have performed a detailed calculation of the branching ratios for the radiative decay Xi v l 
for the lightest neutralino in R-parity violating SUSY models. It is seen that the branching ratio 
can have a maximum value of about 5-10% when R-parity violation has its origin in the so-called 
'bilinear' terms in the superpotential, and is insignificantly small for that induced by 'trilinear' 
terms. This, we point out, can be an interesting way of obtaining characteristic signals for the 
former type of scenario. Our conclusion is that the chances of such radiative decays faking the two- 
photon signals for gauge-mediated SUSY breaking are small, except where the lightest neutralino 
in the latter is close to the kinematic limit of production. It is also seen that the accessible 
region corresponding to the MSSM parameter space does not change appreciably upon subjecting 
the theory to constraints from neutrino mass patterns required by the Superkamiokande data on 
atmospheric neutrinos. 
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Appendix 

Here we outline the actual forms of the various quantities used in the loop calculations in section 
3. 

Expressions for the mixing parameters Al, Bl, Ar, Br for Vl^ -loops in the bilinear R-parity 
violating scenario: 

A L = -^ z n y k2 ~ Z i3 V kl (A.l) 

B L = -^Z n V k2 - Z j3 V kl (A.2) 

Ar = -—j^Z i2 U k2 - Z i3 U k \ + -j=Z i5 Uk3 (A. 3) 

Br = -^Z j2 Uk2 - Z j3 Uki + -^=Zj 5 Uk3 (A. 4) 

where we have j = 2 and i = 1. Similarly for charged scalar loops we have: 



Al = m' (Z i5 V k3 C 2m — Z i2 V k3 C 3m ) + 2Q T tan OwZi^/ k3 C^ m + 

V2{ZnV kl + -^=(Z iA V k2 t<m9 w + Z i3 V k3 )}C lm (A.5) 



A R = -m! Z i2 U k3 C 4m + (Z i3 - Z i4 tan 6 W (1 + 2Q T ))U k3 C 3m + 

V2{Z l2 U kl - -^=(Z l4 U k2 tan 6 W + Z i3 U k3 )}C 2m (A.6) 



Bl = —m Zj 2 U k3 C4 m + (Zj 3 — Zj^tanOwi^ + 2Q T ))U k3 C 3m + 

V2{Z j2 U kl - -^=(Z j4 U k2 tan6 w + Z j3 U k3 )}C 2m (A.7) 



Br = rn'(Zj 5 V k3 C 2m — Zj 2 V k3 C 3m ) + 2Q T tan9wZj4 : V k3 C4 m + 

V2{Z n V kl + —(Zj4V k2 tan W + Z j3 V k3 )}C lm (A.8) 

m' = ^ (A.9) 

(mfo - ^-)2cos/3 

Z's are the neutralino mixing elments, U's and V's are the chargino mixing elements and C's are 
the charged scalar mixing elements. 



where 



Special cases for stop and sbottom: we have considered the effect of mixing between ti and Ir 
as well as between 6l and Ir 
(i) For t-ti triangle: 

A L = (Zf + 2Q t Z u tim9 w )(-sm9 t ) + -— -m t Z xl cos9 t (A.10) 

rriw sm p 

Ar = Z n (-sin^) - 2Q t tan^Z 14 cos^ (A.ll) 

mw sm p 

B L = — ^— Z 21 (-sm9 t ) - 2Q t tan 9 W Z 24 cos 9 t (A.12) 
sin p 

5k = (Z 2 ™ + 2Q t Z 24 tan^)(-sin^) + ^—m t Z 21 cos 9 t (A.13) 

«Hv sm p 

In our convention, t\ is the lightest physical stop and couplings or t 2 can be obtained by replacing 
— sin 9t — ► cos 6>i and cos t — > sin t 
(ii) For 6-61 triangle: 

A L = (-Z7 + 2Q b Z 14 tan6» w O(-sin0 6 ) + 1 -—-m b Z 12 cos 9 b (A.14) 

mvy sm p 

= Tk—^z 12 (-fane b ) -2Q b tan 6 W Z U cos 9 b (A.15) 

ttiw sm p 

-Bl = : — -Z 2 2(-sm9 b ) - 2Q b tan9 w Z 2 iCOs9 b (A. 16) 

mw sm p 

B R = (-Z 2 _ + 2Q b Z 24 tan6» H /)(-sin6> b ) + l — —m b Z 22 cos 9 b (A.17) 

miv sm p 

with 

Zf = Z13 - Z M tan W (A. 18) 



Z 2 = Z23 — Z 2 4tan9w (A. 19) 

In our convention, b\ is the lightest physical sbottom and couplings or b 2 can be obtained by 
replacing — sin 9 b — > cos 0?, and cos 0;, — ► sin #6 

(iii) Degenerate sfermion loops: For up-type squarks we can use the same expression as that for 
the stop loop without the mixing angle corresponding to ti, - tR mixing. Similarly for down- type 
squark or charged sleptons, one can use the expressions corresponding to the sbottom loop modulo 
the mixing factor. Also the charges and the color factors should be properly included. 
In the case of A- type couplings the expressions become (for the down- type squarks): 



A L = O(A') (A.20) 



A R = A(0) 



(A.21) 



B L = ^-5^12 - 2Q d Z u tan 6 W (A.22) 

my/ cos p 

Br = ^—5^12 + 2Q d Z 14 tan 6 W - Zf (A.23) 

mw cos p 

where appropriate generational indices should be considered for the A- and A'-type couplings. 
In the case of charged sleptons one should use the appropriate masses and charge. One must 
remember that in the case of trilinear R-parity violating scenario we have diagonalised a 4 x 4 
neutralino mass matrix. 

The integrals corresponding to the three-point functions are 

If 1 

°0 = — 2 / d <I Tj>. , __.2w_ 35 I Z^777 ,. ^2 I „2n ( A - 24 ) 



'/if — • o 



V 


+ m })(q - 


-P) 2 


+ m 2 s )((q - 


k 2 ) 2 


+ m 2 ) 


5 fa 2 


+ m)){q - 


-P) 2 


+ m 2 )((q- 


k 2 ) 2 


+ m 2 f) 














'v 


+ m 2 )(q 


-P? 


+ m 2 )((q- 


-k 2 ) 2 


+ m 2 ) 



and 



(A.26) 



They are expressed in terms of various form-factors as follows 

C„ = -VyPxi + k^C 12 (A.27) 
Cfxv = P^PuC 21 + k ltl kxvC 22 - {p^k lu + Pvh^Cw + 5^ y C 2i (A.28) 



kip =P»- k 2ph (A.29) 

I 1 = Cia - Cii (A.30) 

I 2 = Co (A.31) 

J 3 = J 4 + J' 4 (A.32) 

/ 4 = C 23 - C 22 (A.33) 
and I' 4 can be obtained from I 4 by m s m/ interchange. For more details on the three-point 



form-factors, see, ref. [28] 



J is analogous to I but the denominator in Eqn.[A.24] will be changed to 

(q 2 + m}){(q - h) 2 + m 2 }{(q - p) 2 + m 2 } 



The integrals which appear in the calculations of M.2 are identical in form to those mentioned 
above but one should replace m s with my/ and m/ with Here m s and m/ are the masses of 
the scalar and the fermion appeared in the loop respectively and is the mass of the chargino. 
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Figure 1: One-loop contributions to the decay Xi — ► v l- I n addition, for every diagram there is a 
counter part with the internal fermion line(s) reversed. 
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Figure 2: Branching ratio for the decay Xi ~^ v l (with SK). The remaining super symmetric 
parameters are chosen as: \i = 200 GeV, B\ = —200 GeV, mj = m^, = 200 GeV, rriq = 600 GeV. 




Figure 3: Branching ratio for the decay Xi ~^ v l (with SK). The remaining super symmetric 
parameters are chosen as: \i = —200 GeV, B\ = 200 GeV, m ; ~ = = 200 GeV, rriq = 600 GeV. 
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Figure 4: Branching ratio for the decay Xi ~^ v l (with SK). The remaining super symmetric 
parameters are chosen as: \i = 500 GeV, B\ = —200 GeV, m ; ~ = rriy = 200 GeV, rriq = 600 GeV. 
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Figure 5: Branching ratio for the decay Xi ~~ * v l- The remaining supersymmetric parameters are 
chosen as: \i = -500 GeV, e = 10 GeV, v 3 = 1 GeV, B x = 200 GeV, m, = mc = 200 GeV, 
nig = 600 GeV. 
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Figure 6: Branching ratio for the decay Xi ~> v l- The remaining supersymmetric parameters are 
chosen as: fi = 500 GeV, e = 10 GeV, v 3 = 1 GeV, B x = -200 GeV, m l = mc = 200 GeV, 
nig = 600 GeV. 
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Figure 7: Branching ratio for the decay Xi v l (with SK), where M\ and M2 are free parameters. 
The value of M2 is taken to be 300 GeV. The remaining supersymmetric parameters are chosen as: 
H = 200 GeV, B x = -200 GeV, m l = mc = 200 GeV, nig = 600 GeV. 
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Figure 8: Branching ratio for the decay Xi ~^ v l (with SK), where M\ and M 2 are free parameters. 
The value of M 2 is taken to be 150 GeV. The remaining supersymmetric parameters are chosen as: 
fi = 500 GeV, B x = -200 GeV, mj = m = 200 GeV, nig = 600 GeV. 
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Figure 9: Branching ratio for the decay Xi ~^ v li where M\ and M 2 are free parameters. The value 
of M 2 is taken to be 300 GeV. The remaining supersymmetric parameters are chosen as: jj, = —500 
GeV, e = 10 GeV, v 3 = 1 GeV, B 1 = 200 GeV, m ! = m = 200 GeV, m,- = 600 GeV. 
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Figure 10: Branching ratio for the decay Xi ~^ v l i n t ne scenario where only A and A' terms are 
present. The remaining supersymmetric parameters are chosen as: \i = —500 GeV, mj = 200 GeV, 
m - = 600 GeV. 



